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Abstract 

Non-attractor models of inflation are characterized by the super-horizon evolution of 
curvature perturbations, introducing a violation of the non-Gaussian consistency relation 
between the bispectrum’s squeezed limit and the power spectrum’s spectral index. In this 
work we show that the bispectrum’s squeezed limit of non-attractor models continues to 
respect a relation dictated by the evolution of the background. We show how to derive 
this relation using only symmetry arguments, without ever needing to solve the equations 
of motion for the perturbations. 


1 Introduction 


The measurement of departures from a purely Gaussian distribution of primordial curvature 
perturbations would give us access to exquisite details about the physics underlying cosmic 
inflation [1—3]. Indeed, different models of inflation predict distinctive deviations from Gaus- 
sianity that are sensitive to the perturbations’ self-interactions, as well as their interactions 
with other degrees of freedom that might have existed during inflation [4,5]. This realization 
has motivated an extensive amount of research over the past decade on the subject of primor¬ 
dial non-Gaussianity, both theoretically [6-13] and observationally [14-17]. Although current 
constraints on non-Gaussianity coming from Cosmic Microwave Background (CMB) observa¬ 
tions remain poor [18-20], future Large Scale Structure surveys [21-23] and Lyman-a Forest 
observations [24-26] promise to substantially improve our knowledge about non-Gaussianity, 
giving us better insights into the elusive nature of inflation. 

Because a complete characterization of non-Gaussianity is somewhat out of reach, it is cus¬ 
tomary to parametrize non-Gaussian departures by defining the bispectrum i? 7 £(ki, k 2 , k 3 ), 
which determines the 3-point correlation function of adiabatic curvature perturbations, hereby 
denoted by TZ, in momentum space 

(^(k!)^(k 2 )^(k 3 )) = (27r) 3 5 3 (k 1 +k 2 + k 3 )5 7 e(k 1 ,k 2 ,k 3 ), (1.1) 

where the Dirac-delta function appears as a consequence of the homogeneity and isotropy of 
the inflationary background. At first order, the functional dependence of B-ji on the three 
momenta ki, k 2 and k 3 is determined by the non-linear evolution of 1Z, parametrized by tree- 
level cubic interactions appearing in the Lagrangian describing its perturbative dynamics. 
As a consequence, in the simplest class of inflationary models —namely single field slow-roll 
inflation— one predicts a scale invariant bispectrum with an amplitude suppressed by the 
slow-roll parameters [9,10], but with a shape determined by the configuration of the momenta, 
restricted to the shell ki + k 2 + k 3 = 0. More exotic single-field (non-canonical) models of 
inflation may predict substantially larger departures from Gaussianity [27,28], some of them 
even showing strong departures from scale invariance [29-33]. These models have been syste¬ 
matically analyzed with the help of the effective field theory of inflation approach [34], which 
allows one to study both canonical and non-canonical models of inflation within a general 
framework without the need of specifying the physics underlying inflation [35-40]. In this 
framework, deviations from canonical inflation are parametrized by the sound speed at which 
curvature perturbations propagate (among other quantities). It is now well understood that 
a suppressed sound speed increases the strength of curvature perturbations’ self-interactions, 
therefore enhancing the amount of primordial non-Gaussianity generated during inflation. 
Models characterized for allowing a suppressed value of the sound speed, with distinctive 
non-Gaussian shapes, include P(A)-inflation [27,28,41], DBI-inflation [42,43], and the low 
energy limit of multi-field inflation with heavy fields [44-46], just to mention a few. On the 
other hand, it is also possible to have large non-Gaussianity appearing as a consequence of 
nontrivial initial states [47-50]. 
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One of the most outstanding discoveries in the study of non-Gaussianity is the so-called 
consistency relation, first reported by Maldacena in ref. [10]. This relation tells us that in the 
squeezed limit (that is, the configuration where the size of one of the three momenta ki, k 2 
and k ;3 is much smaller than the remaining two) the bispectrum is determined by the power 
spectrum Pji and its spectral index n-n in the following specific way 

lim S^(ki,k 2 ,k 3 ) = (n^(/c 3 ) - l)P n {ki)Pn{k^) : (1.2) 

k 3 —>0 

where k t = |k|, and where Pn(k ) and are defined via: 


<^(k 1 )^(k 2 )) ss (k 3 + k 2 ) P n {h), (1.3) 

nn{k) = 1 + ^ In [k 3 P n (k)\ . (1.4) 

Because both the amplitude of the power spectrum P-n and its spectral index are fairly 
well measured [51,52], this consistency relation is regarded as one of the most powerful tools 
available to falsify a large variety of single field inflationary models. It is not even restricted 
to slow-roll models of inflation [53]. Its elegant derivation, which will be reviewed in Sec¬ 
tion 2, combines simple symmetry arguments with simple statistical considerations that are 
independent of the dynamical details of inflation [54-56]. The basic assumption underlying 
its derivation is that, during inflation, the value of every background quantity (for example 
the Hubble parameter H ) is uniquely determined by a single time-dependent parameter (for 
instance, the value of the inflaton vacuum expectation value, were we interested in the case 
of single field inflation). Models respecting this property are called attractor models, and 
they are further characterized by the fact that curvature perturbations freeze after horizon 
crossing. For this reason, a violation of the consistency relation, confirmed by observations, 
would automatically rule out every model of inflation in which the amplitude of curvature 
perturbations remained constant after horizon crossing, encompassing essentially all sensible 
single-field models of inflation. 

One of the preferred ways to constrain non-Gaussianity with the help of CMB observations 
is by defining the so called local /NL-P aram eter, which is related to the squeezed limit of the 
bispectrum through the relation 


f = _A l™ £ft(ki,k 2 ,k 3 ) 
/NL ~ 12k 3 ^o P n (ki)Pn(k 3 )' 


from where one reads /nl — —5 ( n-jz — 1)/12 ~ 0.02 after using the latest constraints on the 
spectral index n-jz. However, projection effects in the measurement of the CMB preclude us 
from directly measuring this predicted value of /nl> implying that the consistency relation is 
in fact equivalent to = 0 [57], which may be compared to the most recent constraints on 
local non-Gaussianity [20], given by = 0.8 ± 5.0. Thus, current observations are rather 
weak in providing a useful assessment of the validity of the consistency relation. Neverthe¬ 
less, the prospects of measuring violations to the consistency relation (1.2) has motivated a 
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fairly big amount of research devoted to the study of models of inflation where curvature per¬ 
turbations are forced to evolve outside the horizon. 1 The most prominent example of such 
models is multi-field inflation, where curvature perturbations have the chance to interact 
with other light degrees of freedom even after horizon crossing, allowing for large deviations 
of the consistency relation [63-65]. Another example is the case of warm inflation [66], in 
which the inflaton remains coupled to a thermal bath. However, it was recently realized that 
curvature perturbations may also grow after horizon crossing in purely single field models 
characterized by a “non-attractor” evolution of their background. For instance, the authors 
of ref. [67] (see also [68]) studied the generation of non-Gaussianity in a class of models known 
as ultra slow-roll inflation [69,70], where the inflaton potential is exactly flat, precluding the 
existence of an attractor regime. In such a model, curvature perturbations grow outside the 
horizon at a dramatic rate, implying a squeezed limit for the bispectrum of the form 

lim -Bft(ki,k 2 ,k 3 ) = -6 Pn(ki)Pn(h), (1-6) 

k 3 -s>0 

corresponding to /nl = 5/2. Given that in these models the power spectrum is almost scale 
invariant n-ji ~ 1, we see that eq. (1.6) represents a flagrant violation of the consistency 
relation depicted in eq. (1.2). Furthermore, it has been argued that large violations of the 
consistency relation may constitute a generic feature of non-attractor models [72,73] which 
may provide an alternative paradigm —to that offered by slow-roll inflation— in order to 
explain the generation and evolution of primordial curvature perturbations. 

After accepting the fact that in this class of models the consistency relation is violated due 
to the super-horizon evolution of curvature perturbations, we ought to ask whether there are 
alternative explanations behind eq. (1.6) other than a brute force computation based on the 
dynamics of the perturbations. The purpose of this article is to clarify the underlying nature 
of the violation of the consistency relation in models with super-horizon evolution. To this 
extent, we will focus our discussion on single field models with non-attractor backgrounds. 
Our aim is to show that models displaying a violation of the consistency relation are still 
restricted to respect well defined relations between the squeezed limit of the bispectrum and 
the power spectrum, similar to Maldacena’s consistency relation. We will show that, in fact, 
the same arguments leading to Maldacena’s consistency relation, if phrased correctly, will still 
reproduce the result expressed in eq. (1.6) found in ref. [67]. In other words, the squeezed 
limit may be determined purely from symmetry and statistical arguments, without the need 
of understanding the details about the dynamics of the specific system under interest. We 
will do this first for ultra slow-roll inflation, and then extend our results to more general non¬ 
attractor models. In this regard, we deduce a general expression determining the bispectrum’s 
squeezed limit, given by 

lim Bft(ki,k 2 ,k 3 ) = -|(4 + r/)P^(/ci)Ffe(£:3), (1-7) 

k 3 —md c% 

1 It is also possible to violate the consistency relation with a non-trivial initial state for curvature pertur¬ 
bations. See for instance refs. [58-62]. 
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where 7] = e/He (with e = —H/H 2 and H being the Hubble expansion rate during inflation) 
and c s is the speed of sound of curvature perturbations. In this expression, the value of 
both and c s will depend on the specific model allowing for non-attractor solutions of the 
background. 

We have organized this article as follows. In Section 2 we review the derivation of the 
consistency relation (1.2) closely following the discussions found in refs. [54,55]. Then, in 
Section 3 we briefly review the model of ultra slow-roll inflation, which is the simplest model 
admitting a non-attractor behavior of the background. Section 4 is devoted to the study of 
the superhorizon behaviour of the inflaton and metric (scalar) perturbations in this model. In 
Section 5 we show how the same arguments used to derive Maldacena’s consistency relation 
can help to set up a modified relation valid for ultra slow-roll inflation, yielding the same 
result derived in ref. [67] by a direct computation. Then, in Section 6 we generalize the 
arguments developed in the previous sections to more general non-attractor models, based 
on P(A)-models of inflation. Finally, in Section 7, we provide our concluding remarks. 

Before commencing, a quick word about units and notation: We shall use natural units 
whereby c = 1, h = 1 and Afpj = 1/87 tGn = 1. In addition, in this paper we formally 
denote the gauge invariant co-moving curvature perturbation by 1Z, which in co-moving 
gauge coincides with the spatial metric perturbation (usually denoted by ip). Furthermore, 
there is a sign difference between our convention for 1Z (and if)) and the convention used in, 
among others, references [10,54,67]. See for example our metric in (2.1). At every stage 
where we cite results obtained in these references, we have accounted for that sign difference. 

2 Review of the consistency relation 

In this section we offer a review of the derivation of the consistency relation (1.2), closely 
following the discussions of refs. [54] and [55] (see also ref. [74,75]). Let us start by recalling 
that the perturbed Friedman-Robertson-Walker (FRW) metric in co-moving gauge may be 
written with the help of the Arnowitt-Deser-Misner formalism [76] as 

ds 2 = -N 2 dt 2 + a 2 {t)e~ 2n 8 ij {dx i + N i dt){dx j + N j dt), (2.1) 

where a(t ) is the scale factor parametrizing the expansion of spatial foliations during inflation, 
and 1Z represents the adiabatic curvature perturbation. In addition, N and N' 1 are the usual 
lapse and shift functions respectively, to be determined in terms of other quantities by solving 
constraint equations. It is customary to define 5N through 

N = 1 + SN, (2.2) 

in which case the background FRW metric is recovered by setting 1Z = 0, SN = 0 and 
N l = 0. Let us notice that an immediate consequence of the metric (2.1) is a symmetry 
under simultaneous rescalements of the scale factor and the co-moving coordinates: 

apt ) —»• a'(t ) = a[t)e AC , dx —»• dx' = dxe~ AC . (2.3) 
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This rescalement does not affect other observable background quantities, such as the Hubble 
expansion rate H = a/a. In addition, notice that the rescalement of the scale factor may be 
absorbed in the curvature perturbation 


TZ -> n' = tz - AC, 


(2.4) 


implying that the equations of motion for TZ must admit at least one solution of the form 
TZ = constant, which reminds us that, after all, we can only measure gradients of TZ. 


2.1 Attractor backgrounds and long wavelength modes 

Our main focus in this section are models of inflation characterized by a background with 
an attractor behavior. These are models where every background quantity is uniquely de¬ 
termined by a single parameter, which may be used as a replacement of time. The simplest 
example of such backgrounds is offered by single field slow-roll inflation, where the back¬ 
ground quickly asymptotes to an attractor trajectory respecting equations of motion of the 
form 

OTA 

3 H 2 = V (<j>), (2.5) 


■ dV 

3 ^ + - = 0 , 


where H = a/a is the Hubble expansion rate. These equations tell us that both H and cf> are 
completely determined by the value of irrespective of the initial conditions. Universes fol¬ 
lowing this feature have been dubbed single-clock universes by the authors of ref. [54], There 
is a direct consequence on the dynamics of perturbations springing out from this behavior. 
First, notice that since we are working in co-moving gauge (that is, scalar perturbations are 
only present in the metric via eq. (2.1)), then every patch of the universe is determined by 
the same value of </>, implying that every patch is characterized by the same Hubble parame¬ 
ter H. Second, recall that the wavelength of any perturbation in a FRW background grows 
proportionally to the scale factor a(t). This means that we can split both TZ and 5N in short 
and long wavelength contributions 


TZ — TZ S T TZc 


5N = 5N S + 5N e , 


( 2 . 6 ) 


where TZp and SNp contain contributions with frequencies smaller than H. A physical observer 
who has access only to short wavelength perturbations will not be able to distinguish the long 
wavelength contributions TZp and 5Np from the background. To understand the consequence 
of this, let us consider a reference time to and insert back the long wavelength contributions 
of the splitting (2.6) in the metric (2.1). We find 

d-s 2 1 = —dt 2 B + a 2 fi(tB,'x)8ijdx l dx : i , (2.7) 


where ts and a e g are given by: 


ts(t,x) = t + St(t,x), 
a e fi(t B , x) = a(t(t jB ))e _7e dhx) ) 


5t(t,x.) = / dt'8Np(t' ,x), 
Jt 0 


( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 
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where t(ts) = t,B — 5t (notice that we have conveniently set St = 0 for t. = f 0 )- Now, we 
wish to examine how the long wavelength fields affect the local background about x at the 
vicinity of t = to. First, notice that a e g- has to be a solution of the same equations of motions 
respected by the original solution a(t). Given that proper time at the x-patch is given by t B , 
the only possibility for a e g to satisfy these equations is that it consists of a evaluated at t B , 
up to a multiplicative constant allowed by the symmetry (2.3) of the background: 

a e ff = a{tB)e ~ AC . (2-11) 


By equating this with (2.9) we find 

a(t)e~ ni = a(t + St)e~ AC . (2-12) 

Then, expanding the right hand side about t , we obtain 

a(t)e-** ~ a{t) [1 + HSt(t, x)] e~ AC ~ a(t)e HSt ^~ AC . (2.13) 

Comparing this result with (2.11) we see that TZ( = AC — H5t, from where we derive that 
1Z( = —H5t — HSNp. Then, by disregarding H5t — HNp(t — to), which is sub-leading 
compared to the other two terms, we conclude that: 

S N ( = -^. (2.14) 

Notice that this is the usual relation obtained by solving the constraint equation for the shift 
5N at linear order. While this result is valid for all wavelengths, our derivation is only valid 
for long wavelength modes. 

Up to this point our arguments have been rather general, and valid for both attractor and 
non-attractor backgrounds. To see the implications of dealing with an attractor background, 
let us go back to eq. (2.10). Recall that, since we are in co-moving gauge, in an attractor 
background there is only one value of (j) characterizing the background at a given time t. In 
particular, the value of the background field at the local patch centered at x must coincide 
with (/>(ts): 

<f>efi(t B ,x) = <f>(tB)- (2.15) 

By comparing this expression with (2.10) we see that 5N^ = 0. In addition, it also implies 
that 

TZf = constant, (2-16) 

is the only allowed behavior for long wavelength curvature perturbations in attractor back¬ 
grounds. In other words, as their wavelengths are stretched, curvature perturbations must 
freeze (implying that the second long wavelength mode must decay). 
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2.2 The consistency relation 


Having established some properties of attractor backgrounds, let us now proceed to derive 
the consistency relation. We want to compute (77(ki)77(k2)77(k3)) in the squeezed limit 
where Aq <C k±,k 2 - In this limit, the long mode k% has left the horizon much earlier than 
those parametrized by k± and & 2 - We can then compute the correlation function between 
77(ki) and 77(k 2 ) in a background renormalized by 77(k3). We begin in position space 
(parametrized by co-moving coordinates) by computing the two-point correlation function 
(7777) (xi, X 2 ) = (77(xi)77(x 2)) in a patch of the universe centered at (xi+X 2)/2 and inflating 
according to the effective scale factor: 

a eS = a(t)e~^, 77; = 77,(^±^). (2.17) 

Now, we may choose to rescale a e g- back to a as long as we properly rescale the co-moving 
coordinates. In other words, we may write 

(7777)b(xi,x 2 ) = (7777) 0 (e-^x 1 ,e-^x 2 ), (2.18) 

where (7777)o denotes the computation of the two-point correlation function in a background 
with an expansion dictated by a(t) alone. However, since (7777)o(xi,X 2 ) is just the two-point 
correlation function computed in the true homogenous and isotropic background, the result 
must depend on the difference |xi — X 2 I. This further implies that: 

(77.77.)s(xi, x 2 ) = (7777) 0 (e'^|xi -x 2 |). (2.19) 

We may now Taylor expand about 77; = 0. Keeping the first two terms of the expansion, it 
is straightforward to find: 

(7777)b(x 1 ,x 2 ) = (7777) 0 (|x 1 - x 2 |) - 77;( Xl + X2 ) --(7777) 0 (|x 1 - x 2 |) H-. 

V 2 / a in |xi — X 2 I 

( 2 . 20 ) 

Next, we can move the expressions to Fourier space. By performing a Fourier transformation 
with respect to xi and X 2 , it is direct to find 

(7777)5(^^2)- (7777) 0 (A: a )+77 < (k<)i^ Tr fei^(Ml > ( 2 - 21 ) 

Kg Cl ill rvg L 

where k^ = ki + k 2 and k s = (ki — k2)/2. As a last step, we may correlate the result 
of eq. (2.21) with 77(k3). By doing this, the first term at the right hand side of eq. (2.21) 
averages out, returning: 

((7777)s(ki, k 2 )77(k 3 )) ~ (77 £ (k^)77(k 3 ))p^ F [k 3 s Pn(k s )]. (2.22) 

The left hand side of this expression gives us back (77(ki)77(k2)77(k3))k 3 ->.0) finally leading 
to the desired result 

lirn B n ( ki,k 2 ,k 3 ) = {n n {k 3 ) - Pnih)Pnih), (2.23) 

k 3 —>0 
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where the spectral index is identified with: 


nn{k) = 1 + In [k 3 Pn(k)] . (2.24) 

The consistency relation tells us that any measurement of /nl larger than 1 — ~ 0.04 will 

rule out every model in which the superhorizon perturbations R£ freeze, i.e. every attractor 
adiabatic single field slow-roll inflation model. On the other hand, in the next section we will 
see that in the non-attractor model of ultra slow-roll inflation, the modes 7Z( evolve rapidly 
outside the horizon. Then it does not come as a surprise that the consistency relation breaks 
down, as reported in [67]. 

3 Non-attractor backgrounds: Ultra slow-roll inflation 

In this section we examine some aspects of non-attractor models, where the evolution of the 
background keeps some memory of the initial conditions. To keep our discussion simple we 
will focus our attention on the particular case of ultra slow-roll inflation studied in refs. [69, 
70]. We shall examine more general non-attractor backgrounds in Section 6. 

3.1 Ultra slow-roll inflation 

Ultra slow-roll inflation was first introduced in ref. [69] and further worked out in ref. [70]. In 
its simplest version, ultra slow-roll inflation is realized by an exactly flat scalar field potential 
of the form 

V(<fi) = Vq, (3.1) 

where Vo is a constant. Evidently, this potential does not constitute a very realistic choice (see 
for example the discussion in [68]) but it allows us to study the evolution of perturbations in 
backgrounds that are dramatically different from those encountered in conventional slow-roll 
inflation. Let us notice that eq. (3.1) automatically implies that the theory is invariant under 
the shift symmetry 0 —> 0' = 0 + Ac j> in addition to the rescaling of the scale factor examined 
in Section 2. We will come back to the consequences of this symmetry in a moment. The 
background equations of motion of the system are 

0 + 3770 = 0, (3.2) 

QH 2 = <j ) 2 + 2U 0 . (3.3) 

The first equation already tells us that there is no such thing as an attractor behavior 
dominated by the friction term 3770 except for the trivial solution 0 = 0. By adopting the 
notation H = 77(0), these equations can be combined into a single equation given by: 

377 2 - 2(77 / ) 2 = Vq. (3.4) 

The solution to this second equation gives us 77 in terms of 0 once we have adopted a set 
of boundary conditions. For definiteness, we choose these conditions in such a way that 



3 H 2 = Vq at 4> = 0. This condition implies: 


H(#) = 



(3.5) 


To continue, we may choose to study inflation in the range <f> < 0. Then, the equation of 
motion for (p(t) may be integrated once to give: 


<t> 



(3.6) 


These equations permit us to compute various background quantities in terms of (f>. For 
instance, the scale factor a{4>) is found to be 


a{4>) = a 0 



(3.7) 


whereas the slow-roll parameters e = —H/H 2 , = e/(eH) and £ = r]/(r)H) are respectively 

found to be given by 


e((j>) = 3tanh 2 ^ \/3/2 (j^j , 

(3.8) 

r}{4>) = —6 cosh -2 ^-y/3/2 (frj , 

(3.9) 

£(</>) = 6tanh 2 (y/3/2 (f>\ . 

(3.10) 


These results, in combination with eq. (3.6), show that in the limit t —>• +oo the background 
parameters have the following asymptotic behavior: 


e —>■ a 6 , rj —>• —6, £ —>• a 6 . 


(3.11) 


Thus we see that this inflationary background is characterized by a substantially large value 
of r), offering a large departure from conventional slow-roll inflationary models. In what 
follows, we will restrict our analysis to the particular case of e <C 1, which is necessary to 
reproduce a scale invariant power spectrum (See appendix A). 


3.2 Long wavelength perturbations 

Just as we did in Section 2.1, we may infer the behavior of long wavelength modes in co- 
rnoving gauge by examining the symmetries of the non-attractor background at hand. To 
start with, recall that the ultra slow-roll background of eq. (3.1) is characterized by the 
symmetry: 

(f) —>• ft = 0 + A cj). (3-12) 

Let us emphasize that, as opposed to the case of attractor backgrounds, in ultra slow-roll 
inflation there are no solutions uniquely linking (background quantities like) H and cj). A 
variation of the initial conditions leads to a variation of the relation between H and <j). This 
implies that, as the wavelength of perturbations are stretched by the expansion of space, the 
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long wavelength contribution to of eq. (2.6) may modify the background value of 4> felt 
by sub-horizon modes, just as in eq (2.10), in a way consistent with the symmetry (3.12). In 
other words, this time, in addition to eq. (2.11), we may have 

= 4>(t B ) + A</>, (3.13) 


where A<j) is a constant incorporating the effects of long-wavelength perturbations via the 
shift function Ng (but should not be confused with perturbations of the field </>, which in 
co-moving gauge are turned off). From (2.10) we see that the previous equation is equivalent 
to 

4>{t) = (j)(t + St) + A (f>. (3-14) 


This implies that to linear order: 


St = 


Acp 

4 > 


Then, using this result in combination with (2.11) and (2.9) we obtain 


(3.15) 


Tie 


AC + ^-Acf>. 

$ 


(3.16) 


Recalling eqs. (3.5) and (3.6) of the previous section, we see that the asymptotic behavior of 
1Z{ soon or later will be dominated by A<f> which implies a growing mode of the form: 


7Zf> ~ a 3 . 


(3.17) 


That is, ultra slow-roll inevitably contains super-horizon evolution of curvature perturbations. 
We are for sure not the first ones to achieve this result, see for example [67,70]. However, 
we would like to emphasize that in this case it has been exclusively deduced with help of 
symmetry considerations. 


3.3 A violation of the consistency relation? 

The result expressed in eq. (3.17) tells us that the argument used to derive the consistency 
relation involving the rescalement of the co-moving coordinates cannot be repeated without 
carefully taking into account the additional contribution from (j) e g felt by sub-horizon modes. 
Indeed, a brute force computation of the squeezed limit of non-Gaussianity leads to a violation 
of the consistency relation [67] in the form: 

lim H^(ki,k 2 ,k 3 ) = —6P^(k 1 )P 7 ^(k 3 ). (3.18) 

k 3 —>o 

This result signals a violation of the consistency relation (2.23) which in slow-roll attractor 
backgrounds prescribes 


lim Rft(ki,k 2 ,k 3 ) 

k 3—>0 


-(?? + 2e)P^(ki)P^(k 3 ), 


(3.19) 
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where we have used the standard relation between the spectral index and the slow-roll param¬ 
eters 1 — n-jz = i]+ 2e. In this case, the mere fact that a model of adiabatic single field slow roll 
inflation can produce an order-one non-Gaussianity (one gets /nl = 5/2) shows the apparent 
breakdown of the consistency condition 2 . Thus, as the consistency relation is supposed to be 
valid for all adiabatic single-clock models of inflation, the model of ultra slow-roll inflation 
seems to provide a counterexample. However, in section 5 we will show how a generalized 
version of the consistency relation still holds for the model of ultra slow-roll inflation. 

4 Freezing of superhorizon perturbations 

Before analyzing the squeezed limit in the context of non-attractor models, let us briefly 
analyze the non-linear relation between curvature perturbations and inflaton perturbations 
in the context of ultra slow-roll. Up to this point, we have been working with adiabatic 
perturbations 1Z in co-moving gauge, where inflaton perturbations are absent. More generally, 
TZ is a gauge invariant quantity that reduces to spatial curvature perturbations of the metric 
i/) under the condition that inflaton perturbations 8(f) vanish: 

TZ\sij>=o = 'if’- (4.1) 

Alternatively, we may work in flat gauge, and define a gauge invariant perturbation Q that 
reduces to inflaton perturbations 8<f> under the condition that spatial curvature perturbations 
of the metric if> vanish: 

QU=o = S</>, (4.2) 

It is then possible to find (see for instance [10]) that 1Z and Q are non-linearly related, up to 
quadratic order, in the following way: 

H „ n H 2 ~ H . . . 

K = Z Q + W e -^ ee+ -’ (43) 

or equivalently 

Q = ^n + ^iz 2 + ^izn + ---, (4.4) 

where the ellipses “• • • ” represent terms with spatial gradients. Now, as we have seen, in 
conventional slow-roll inflation 1Z freezes at long wavelengths, implying that TZp —>• 0 fast 
enough to imply the relation: 

Slow-roll: Qi = ±K e + ^1$, n t = jQ t -^Ql (4.5) 

2 Note that the disagreement between (3.18) and (3.19) continues to hold in the decoupling limit: inserting 
e —¥ 0 and r/ —>• —6 still produces a sign difference between the two. It is of course not completely fair to 
compare the two results in such a way, since the derivation of the result (3.19) supposes e, r/ <C 1. 
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However, in the case of ultra slow-roll we have that 7 Zg oc a 3 , from where we read 77/ = 3HTZp, 
giving us back the relation 

Ultra slow-roll: Qp = j-TZp - j77 2 , ^ = T 2 * + (4.6) 

H H 4 i f (j) 1 4 

where we have used i] = —6 + 2e. To exploit these relations, let us examine what happens if 
we were to expand a few quantities encountered in Section 2 up to second order in the long 
wavelength perturbations. First, if we expand (2.12) about t up to second order we find: 

a{t)e~ nt ~ a(t) 1 + HSt(t,x) + ^H 2 ( 1 - e)5t 2 e~ AC ~ a{t)e H5t ^^ eH2st2 ~ AC . (4.7) 

This relation tells us that 

77/, = AC - HSt(t , x) + ^eH 2 5t 2 . (4.8) 

On the other hand, by expanding (3.14) about t to second order, we find 

4>{t) ~ (pit) + (j)5t + -^7 2 + 

~ ^)(7) + 057 — - H<j)8t 2 + A(j), (4.9) 

where we used the background equation of motion <f> + 3H(f> = 0. This equation allows us to 
deduce St up to second order in A<j>: 

* = + (4 ' 10) 

Plugging this result back into eq. (4.8) we finally deduce 

Up = AC+ S-Acp - ^l^Acp 2 , (4.11) 

<p Z (p z 

which, after noticing that r] = —6 + 2e, we see that is precisely consistent with (4.6) once 
we replace Qp —> Acft (and disregard AC). Thus, we see that Q must freeze in order to have 
a long wavelength limit consistent with the symmetries of the background. In fact, a brute 
force computation reassures us that Q indeed freezes [70]. For instance, to second order, the 
action for Q is found to be given by: 

S Q = lf d4xa3 [S' - 2eV oQ 2 ] • (4.12) 

Notice the presence of the mass term proportional to e which in fact breaks the shift symmetry 
tp —i> <f>' = (f) + A<f) shared by the background. However, since t < 1 and e oc a~ b it is still 
compatible with the freezing of Q. To see this, we may perform the following field redefinition 
u = Q/f{4>) with / a function of the background field (f>. The quadratic action in terms of u 
then reads 

Si 2) = \j d 4 xa 3 (f 2 u 2 + <P 2 [//" + • • • ] u 2 ) , (4.13) 
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where the ellipses “• • • ” represents other terms containing / and its derivatives with respect 
to (j>. To find the freezing solutions for u, the part in square brackets above needs to vanish. 
This gives us a second order equation for /. One of the solutions will correspond to u = 1Z 
(there is still a constant mode in 1Z, but it is subdominant compared to the other mode which 
grows as a 3 ). The other solution corresponds to u = Q( 1 + 0(e)), which is just another way 
of stating that (the dominant mode of) Q freezes up to all orders in the perturbations in the 
decoupling limit e —>• 0. 

5 A consistency relation for ultra slow-roll inflation 

Our next goal is to compute (77(ki)77.(k2)77.(k3)) in the squeezed limit k;$ <C k\,k^ for ultra 
slow-roll inflation, employing symmetry arguments as in the case of attractor models. As 
we have discussed in the previous sections, the main challenge is to take into account the 
super-horizon evolution of curvature perturbations 1Z. To proceed, we adopt the flat gauge, 
whereby spatial curvature perturbations ijj are turned off. In this gauge the metric line 
element takes the form 

ds 2 = —N 2 dt 2 + a 2 {t)5ij{dx l + N l dt)(dx J + N^dt), (5-1) 

and scalar perturbations enter as excitations of the inflaton field: 

<Kx,t) = 4>{t) + Q(x,t). (5.2) 

As argued in the previous section, the symmetry of the theory implies that Q(x, t) asymptotes 
to a constant on super-horizon scales. However, given that we derived this result in co-moving 
gauge, it is instructive to show how this result is recovered in flat gauge. This time, we must 
split both Q and N into short and long wavelength contributions of the form 

Q = Qs + Qe, N = N s + Ng. (5-3) 

Then, a patch of the universe centered at x will be characterized by a background with an 
effective field (j) e s and an effective scale factor a e g given by 

(j> eS (t B ,x) = + Qe(t,x), (5.4) 

a eff (f s ,x) = a(t(t B )), (5.5) 

where t(t B ) = t B — St and is given by eq. (2.8). Now, notice that the only form for a e ff 
consistent with the symmetries of the background is a e g = a(ts)e AC ', with A C a constant. 
This implies that: 

a(t) = a(t + St)e~ AC ~ ae~ AC+HSt . (5.6) 

Then, we see that St = AC/H, or equivalently, SNg oc dH~ l /dt. As a consequence, SNg 
vanishes quickly 

SNi oc e -> a -6 , (5-7) 
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and the only long wavelength contribution to (j) e s is via Q? as: 


0 e ff(i,x) = 4>(t) + Qe{t,x). 


(5.8) 


To continue, given that Qi —>■ constant, we may compute the three-point correlation func¬ 
tion (Q(ki)Q(k 2 )Q(k; 3 )) in the squeezed limit k% <C k \, /c 2 using the same arguments of 
Section 2.2. To this extent, we start by computing the two-point correlation function 
(QQ)(xi,X 2 ) = (Q(xi)Q(x 2 )) in a patch of the universe, centered at (xi + x 2 )/2 which 
feels a background given by 

4>eS = 4>(t) + Qe ( Xl ^ X2 ) • ( 5 - 9 ) 

The two-point correlation function (QQ)b(xi, x 2 ) computed in this background may then be 
Taylor expanded about the background without the long wavelength contributions, giving us 
back 


(QQ)b(xI, x 2 ) ~ (QQ)o(|xi 
^ (QQ)o(|xi 


x 2 | ) + a(ii±^ 

x 2 | ) + a(^ 


4-(2S)b(x l,x 2 ) 

L dQe 

-jl(QQ) o(xi,x 2 ), 


(5.10) 


where we have traded the derivative with respect to Qi with a derivative made with respect 
to <f>, the background value of the scalar field. The next steps are straightforward. First, we 
move to Fourier space by performing a transformation with respect to xi and x 2 to obtain 

(QQ)b( ki,k 2 ) ~ (QQ) 0 (k a ) + Qe (k<) ^rP Q (k s ), (5.11) 


where k^ = ki + k 2 and k s = (ki — k 2 )/2. Then we correlate the result with a mode Q(k, 3 ). 
This gives us the squeezed limit: 


lim (Q(k 1 )Q(k 2 )Q(k 3 )) = (2vr) 3 5 3 (k 1 + k 2 + (5.12) 

k 3 —>-0 Ocp 


To obtain a more explicit result we need an expression for the power spectrum Pg(fc) in terms 
of the background quantities. 3 To this extent, it suffice to know that the power spectrum is 
proportional to H 2 (see for example [70]) 


P Q (k) oc H 2 . 


(5.13) 


Then, knowing that H = H((f>), we may write d^PtQ^k) = 2PQ(k)d c f > H/H, which together 
with the fact that d^H/H = \J e/2 (easily derived from eq. (3.5)) leads to 

lim (Q(k 1 )Q(k 2 )Q(k 3 )) = s/2e (27r) 3 5 3 (k! + k 2 + k 3 )P Q (k 1 )P s (k 3 ). (5.14) 

k 3 —>0 

3 Notice that we only need the amplitude of the power spectrum in terms of the background fields evaluated 
at a pivot scale, and we need not to worry about its scale dependence. 
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We can now use this result to find the three-point function of 1Z. In eq. (4.4) we already 
found how to rewrite Q in terms of 1Z and R to quadratic order. Using that in ultra slow-roll 
we have 1Z ~ a 3 while 77 ~ — 6 , then we can write 

2 = 4 In-^U 2 ] +••• . (5.15) 

H l 4 

where the ellipses “• • • ” stand for spatial gradients, which in the present discussion may be 
disregarded. Note that the 1Z term effectively changes the sign of the correction to the first 
order result. This is to be equivalent to the observation made in [67] that neglecting the 
contributions of the decaying mode (i.e. the i?-effects) changes the sign of the final result 
for the computation of /nl- To continue, notice that in general, for a relation of the form 
/ = g + A g 2 (see for example [ 10 ]) one finds: 

(/( x )/(y)/( z )) = (g(x)g(y)g(z)) + 2A [{g(x)g(y)){g(x)g(z)) + cyclic] . (5.16) 

Then, by choosing / = Q/V 2e, g = 7Z and A = — 77 /4, we obtain 

(77(x)77(y)77(z)) = (2e )” 3/2 (Q(x)Q(y)Q(z)) + | [(7Z(x)7Z(y))(7Z(x)7Z(z)) + cyclic] . (5.17) 

To finish, we just need Fourier transform this expression, take the squeezed limit and replace 
the three-point correlation function by our previous result (5.14), to obtain: 

lim (77(k 1 )77(k 2 )77(k 3 )) = (2e + rj) (2 7 r) 3 (5 3 (k 1 + k 2 + k 3 )P^(ki)P^(k 3 ). (5.18) 

k 3 —>0 

This coincides with the result (3.18) obtained in [67] in the decoupling limit e —>• 0 and 
77 —>• — 6 . We wish to stress, once more, that to derive this result we never needed to deduce 
the explicit solution of 7Z nor Q. Instead, we have only used symmetry arguments to relate 
long wavelength modes to the background. It is true, however, that even without explicitly 
deriving it, to deduce our results we have assumed that eventually the faster growing solution 
(the one proportional to a 3 ) will dominate the other one (the constant one). 

The result (5.14) is our proposed modification of Maldacena’s consistency relation, valid 
for the model of ultra slow-roll inflation. However, we have to admit that strictly speaking it 
is not a true consistency relation. The original consistency relation relates three observables: 
power spectrum, bispectrum and spectral index. In this case, the power spectrum and 
bispectrum are related via the factor ( 2 e + 77 ), which is not a direct observable (although for 
ultra slow-roll it simply asymptotes to — 6 ). 

6 More general non-attractor backgrounds 

Before concluding, let us show how our arguments may be implemented in more general non¬ 
attractor backgrounds. To this extent, let us consider inflation in the context of P(Al )-models, 
where the Lagrangian is given by a general function of the inflaton field in the following form 

C = P(X,<!>), X = -±(d<t>) 2 . (6.1) 
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We remind the reader that even in these more general models, we still restrict ourselves to 
models in which e is zero or decaying very rapidly. The background equations of motion are 
found to be given by 


(2XPxx + Px)(f> + 2 XPx<i> + 2>HPx(t> — P^ = 0, 

3 H 2 = 2 XP x - P, 


( 6 . 2 ) 

(6.3) 


where P x = dxP, P^ = d^P, and now X = <fi 2 / 2. 
convenient equation 

H = -XP x , 


One can also deduce the additional 

(6.4) 


which in turn, tells us that 


_ XP x 
e H 2 ' 


(6.5) 


These equations of motion admit both attractor and non-attractor backgrounds. Of course, 
we will be interested in non-attractor backgrounds as long as they satisfy the condition e <C 1 
(in order to have a quasi-de Sitter geometry). Our next challenge is to identify a symmetry 
of the Lagrangian allowing us to deduce how the perturbations evolve after horizon crossing. 
A transformation of the Lagrangian under such an alleged symmetry will satisfy 


A C = P x XX + P 4 A0 = 0, (6.6) 

from where one reads P^ = —P x AX/ Acj). Inserting this result back into the background 
equations of motions we find: 


(2 XP xx + Px)4>A<t> + 2XP x<p Acl> + 3Hj>P x Acf> + P x AX = 0. (6.7) 


By noticing that AX = cj)Acp = 
solution is found to be given by 


(j)d(A(f>)/dt, eq. (6.7) may now be integrated once. The 


A <j) 


AC 

a 3 V2XPx ’ 


( 6 . 8 ) 


where AC denotes an integration constant. Thus, a transformation of the inflaton field of 
the form 4> —>• cj) 1 = 4> + A<f> with Acj) given by (6.8) constitutes a symmetry of the background. 
Notice that in the specific case of ultra slow-roll, examined in Section 3, we have P = X — Vo, 
P x = 1, and \J2X = cj) oc a 3 , and eq. (6.8) is reduced to a shift symmetry Acj) = constant, 
as it should. On the other hand, in the case of a slow-roll background, we recover Acj) ~ a -3 , 
reminding us about the attractor nature of slow-roll inflation. 


6.1 Symmetry and long wavelength perturbations 

To explore the consequences of eq. (6.8) on the perturbations, let us stick to the metric (5.1) 
written in flat gauge. Then, repeating the arguments of Section 3.2, we see that 5N 4 oc e. 
Then, if the non-attractor background is such that e —>• 0 quickly, eq. (6.8) informs us that 
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Ng —> 1, and the value of the effective scalar field (j) e g at a given patch of the universe centered 
at x will be given by: 

ies(t, x) = </>(t) + Qg(t , x). (6.9) 


Then, because this form of the field has to be consistent with the symmetry underlying 
eq. (6.8), we conclude that the long wavelength behavior of the scalar perturbation is given 
by: 


Qe — 


C 


a :i V2XP x 


( 6 . 10 ) 


This in turn tells us that there is a perturbation P, proportional to Q, that freezes outside 
the horizon, and is given by: 

F = a 3 V2XP x Q. (6.11) 


On the other hand, let us recall that, up to second order in the fields, Q is related to the 
adiabatic curvature perturbations 1Z in the following way 


Q 


i. 

H 


n + n 2 + 

4 C 2 


1 

Pci 


nil 


+ • • •, 


( 6 . 12 ) 


where this time we have accounted for the presence of the speed of sound c s of curvature 
perturbations, which for the models examined up to this point had a value c s = 1. The 
explicit value of c s in terms of other background quantities is: 


c 


2 

s 


Px 

Px + 2 XP X x 


(6.13) 


Equation (6.12) further implies that: 


P = 2 P He 


iz + —n 2 h _ — 

/C+ 4 C 2 Hc 2 s 


nn 


+ 


(6.14) 


where we have used the identity of eq. (6.5). Now, at linear order we have n = J-/2a 3 He. 
Therefore, in the long wavelength limit Tig oc 1/2 a 3 He, from where we deduce that at linear 
order 

n e = -H( v + 3)n e , (6.15) 

which is obtained after using the background equations of motion, and the fact that f<l. 
Then, inserting this result back into eq. (6.14) we obtain 


J-g = 2 a 3 He \jlg + A n 2 , 


(6.16) 


where 

A = -£j( 4 + „), (6.17) 

Equation (6.16) is what we need to derive the squeezed limit for the bispectrum of n. We 
examine this derivation in the following discussion. 
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6.2 A general squeezed limit for non-attractor inflation 

In what follows, we derive a general expression for the squeezed limit of non-Gaussianity, 
valid for general non-attractor models. To start with, by using eq. (5.16) with / = J r /2 a 5 He 
and g = 77, we obtain: 

(77(x)77(y)ft(z)) = 8 ffl 9 ff 3 e 3 (^(x)^)^ 2 )) - 2 A[(77(x)77(y))(77(x)77(z)) + cyclic] . 

(6.18) 

Then, we may compute (J r (ki)J r (k 2 )J : "(k 3 )) following the same procedure employed in Sec¬ 
tion 5 for the case of ultra slow-roll. In this case, we obtain 

lim (J r (k 1 )J r (k 2 )J r (k 3 )) = (27r) 3 5 3 (k 1 + k 2 + k 3 )Pr(h) -j^P^h) , (6.19) 

k 3 ->-0 \pJ-l Jo 

where Pjr(k) represents the power spectrum for T computed in a background that is renor¬ 
malized by the long wavelength contributions Jy. This result may be reexpressed in terms 
of the curvature power spectrum as: 

lim (J’(k 1 )J’(k 2 ) J-(k 3 )) = 16 h T a V2 H A e\ 27r) 3 5 3 (k 1 + k 2 + k 3 )Pn(h)Pn(k 3 ). (6.20) 

k 3 -^0 

where we have defined: 

hr=\^lnP^] . ( 6 . 21 ) 

L J o 

Now, putting together eqs. (6.20) and (6.18) in momentum space, we finally arrive to 

lim (7e(k 1 )^(k 2 )^(k 3 )) = 2/^a 3 iA(27r) 3 5 3 (k 1 + k 2 +k 3 )P n (k 1 )P n (h) 

k 3 ->0 

-4A (2vr) 3 <5 3 (k 1 + k 2 + k 3 )P^(fci)P^(fc 3 ). (6.22) 

To simplify this expression, we notice that the first term of the right hand side should vanish 
in the long wavelength limit. Indeed, on the one hand the overall factor 2 a 3 He decreases 
quickly for non-attractor models (otherwise curvature perturbations would not grow). On 
the other hand, because P freezes in the long-wavelength limit, Pjr must tend to a constant, 
independently of the value of Jy (as in the case of ultra slow-roll). These considerations lead 
to our final result valid for general non-attractor models: 

lim ^(ki,k 2 ,k 3 ) = \{A +ri)P n {k 1 )P' R .{k 3 ). (6.23) 

k 3 —>0 Cg 

The overall factor at the right hand side is in general time-dependent, and therefore must 
be evaluated at the end of the non-attractor phase, which could coincide with the end of 
inflation, or could be the beginning of the attractor phase (see refs. [67,71] for a discussion 
on the phenomenological feasibility of this class of models). To finish, let us notice that 
our result reduces to the already known answer found in ref. [71], where a specific family of 
non-attractor solutions were studied. 
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7 Discussion and conclusions 


We have studied the violation of the non-Gaussian consistency relation within single field non¬ 
attractor models of inflation, characterized by the fact that curvature perturbations do not 
freeze after horizon crossing. Our analysis was based purely on symmetry considerations: we 
found that the same arguments leading to the universality of the the non-Gaussian consistency 
relation for attractor models of inflation, may be reproduced to deduce a relation valid for 
non-attractor models. To achieve this, it was important to notice that while curvature 
perturbations 1Z do not freeze after horizon crossing, other perturbations do (Q in the case 
of ultra slow-roll, and J- in the case of more general non-attractor models). As a result, 
it is possible to further understand the violation of the consistency relation as a natural 
consequence of the super-horizon evolution of the modes, dictated by the background of 
the model, independently of the model being an attractor or not. Our results agree with 
those of previous analyses. In particular, we have re-derived the violation to the consistency 
relation in the context of ultra slow-roll inflation reported in ref. [67], and we have found a 
general expression for the squeezed limit of non-Gaussianity for non-attractor models which 
generalizes the results found in ref. [71]. 

Our results offer a rationale to understand potentially large violations to the consistency 
relation, telling us how to relate the size of the violation to the evolution of the inflationary 
background. Given that future large scale structure surveys will offer substantially better 
constraints on non-Gaussianity than those currently available from CMB experiments, it is 
particularly timely to understand the difference between alternative mechanisms to violate 
the consistency relation. In this regard, one particularly interesting task ahead is to charac¬ 
terize the squeezed limit of other n-point correlation functions within non-attractor models. 
Just as in the case of single-field attractor models, it should be possible to derive specific 
relations among different n-point correlation functions, consistent with the symmetries of the 
inflationary background [77-83]. 
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A Evolution of modes in ultra slow-roll inflation 

In this appendix we explicitly compute the behavior of curvature perturbations for both 
standard attractor slow-roll inflation and for non-attractor ultra slow-roll inflation, to linear 
order. These computations are of course complementary to the derivations based on sym- 
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metry arguments found in Sections 2.1 and 3. However, we will use them to re-derive the 
freezing of the modes of Q on super-horizon scales 

Let us start by quickly reviewing the freezing on super-horizon scales of curvature per¬ 
turbations 77. In single-field models it is common lore to define the canonically normalized 
perturbation v in terms of 1Z as: 


v = a—TZ = a —+ 5<f> 

H \ H 


In terms of its Fourier modes v k the equation of motion is 


v k ~ ~ v k + k 2 v k = 0, 


z = 


a<p 

7r 


(A.l) 


(A.2) 


where the prime ' denotes a derivative with respect to conformal time r. To zeroth order 
in the slow-roll parameters, in standard slow-roll inflation we find that z"/'z = 2/r 2 , where 
conformal time may be written as r = 1 /(aH) where H is constant. Then, in the long 
wavelength limit we obtain: 


v k - ~ v k = 0, 

T- 


Vk = ClT 2 


C 2 Cl 

H-oc —^ + C20. 

r a z 


(A.3) 


To zeroth order in the slow roll parameters the factor < j)/H in (A.l) does not change in 
time, so we find that the modes 7 Z k indeed freeze on super-horizon scales. A more detailed 
computation shows that this freezing is exact up to all orders in slow-roll parameters. For 
our discussion the rough sketch above suffices. 

How does this situation change in ultra slow roll-inflation? In this case, e = cj) 2 /(2H 2 ) 
falls down as a -6 and rj asymptotes to —6. Surprisingly, to zeroth order in e and all orders 
in r] (compatible with the decoupling limit), we still find z"/z = 2/r 2 and r = 1 /(aH). To 
be precise, we still have 


z 


u 


z 



1 

4 



(A.4) 


as ij —> —6. Therefore, the solution found in (A.3) for the modes v k is still valid. Things 
change, however, in the conversion to 7 Zk, for which we now find the long wavelength limit: 


7 Zk = 


1 H 

— r v k ~ a 

a cj) 


3 


(A.5) 


As we already found in (3.17) by our symmetry argument, in ultra slow-roll inflation the 
modes 1Z k do not freeze on super-horizon scales. However, from the above it is clear (at least 
up to zeroth order in e and all orders in rj), that the modes of a related perturbation 


Q = 


H 


V’ + 



(A.6) 
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do freeze on super-horizon scales. Note that Q is as perfectly gauge invariant as 1Z, since the 
two only differ by a background quantity. Therefore, both variables share the same scalar 
spectral index n-jz = uq. (The fact that the modes 1Z}. keep evolving after horizon crossing 
does not take away that at any desired moment in time the relative differences in power 
between modes of different wavelengths are encoded in ri-jz ■) 
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